Photodissociation and formation of an ion-pair in CH2 FCl (HCFC-31).

Although CH2 FCl (HCFC-31) recently became of great atmospheric importance, studies concerning its excited states are almost nonexistent. Several excited singlet states were studied (valence nσ* and Rydberg n3s, n3p, σ3s, and σ3p) through highly correlated multireference configuration interaction with singles and doubles, including extensivity correction. Comparison with the states of CH3 Cl indicates a strong influence of the F atom. Potential energy curves suggest formation of an electrostatically bound complex that relaxes to a hydrogen-bonded contact ion-pair (HBCIP) which can decay yielding CH2 F + Cl or to the ground state minimum of CH2 FCl. The HBCIP has a dipole moment of 9.57 D, a CI wavefunction described as 0.65ionic + 0.20biradical and it is strongly bonded by 4.72 eV. Its H bond has characteristics of moderate and strong H bonds. The simulated absorption spectrum confirms the nσ* assignment for the first and suggests the n3s + n3pσ assignment for the second band.

Quantification of the Ionic Character of Multiconfigurational Wave Functions: The Qat Diagnostic.

The complete active space self-consistent field (CASSCF) method is a cornerstone in modern excited-state quantum chemistry providing the starting point for most common multireference computations. However, CASSCF, when used with a minimal active space, can produce significant errors (>2 eV) even for the excitation energies of simple hydrocarbons if the states of interest possess ionic character. After illustrating this problem in some detail, we present a diagnostic for ionic character, denoted as Q at, that is readily computed from the transition density. A set of 11 molecules is considered to study errors in vertical excitation energies. State-averaged CASSCF obtains a mean absolute error (MAE) of 0.87 eV for the 34 singlet states considered. We highlight a strong correlation between the obtained errors and the Q at diagnostic, illustrating its power to predict problematic cases. Conversely, using multireference configuration interaction with single and double excitations and Pople's size extensivity correction (MR-CISD+P), excellent results are obtained with an MAE of 0.11 eV. Furthermore, correlations with the Q at diagnostic disappear. In summary, we hope that the presented diagnostic will facilitate reliable and user-friendly multireference computations on conjugated organic molecules.

Surface hopping modeling of charge and energy transfer in active environments.

An active environment is any atomic or molecular system changing a chromophore's nonadiabatic dynamics compared to the isolated molecule. The action of the environment on the chromophore occurs by changing the potential energy landscape and triggering new energy and charge flows unavailable in the vacuum. Surface hopping is a mixed quantum-classical approach whose extreme flexibility has made it the primary platform for implementing novel methodologies to investigate the nonadiabatic dynamics of a chromophore in active environments. This Perspective paper surveys the latest developments in the field, focusing on charge and energy transfer processes.

Non-Kasha fluorescence of pyrene emerges from a dynamic equilibrium between excited states.

Pyrene fluorescence after a high-energy electronic excitation exhibits a prominent band shoulder not present after excitation at low energies. The standard assignment of this shoulder as a non-Kasha emission from the second-excited state (S2) has been recently questioned. To elucidate this issue, we simulated the fluorescence of pyrene using two different theoretical approaches based on vertical convolution and nonadiabatic dynamics with nuclear ensembles. To conduct the necessary nonadiabatic dynamics simulations with high-lying electronic states and deal with fluorescence timescales of about 100 ns of this large molecule, we developed new computational protocols. The results from both approaches confirm that the band shoulder is, in fact, due to S2 emission. We show that the non-Kasha behavior is a dynamic-equilibrium effect not caused by a metastable S2 minimum. However, it requires considerable vibrational energy, which can only be achieved in collisionless regimes after transitions into highly excited states. This strict condition explains why the S2 emission was not observed in some experiments.

Acid Enriched with Methanol: Formation of a Prebiotic Cluster in the Interstellar Medium.

Organic molecules are a potential source of prebiotic chemistry in the interstellar medium (ISM). Methanol (MetOH) is a very important source of more complex molecules. H3 O+ (aq) and Cl- (aq) are fundamental to living organisms and can be generated in the ISM from the dissociation of HCl with just four water molecules, yielding the (H3 O)+ (H2 O)3 Cl- ion-pair. Here, a detailed mechanism, based on density functional theory (DFT) and ab-initio (2nd order Mϕller-Plesset perturbation theory, MP2) calculations, is suggested for the substitution reactions of these water molecules by MetOH. The time required for formation of an appreciable amount of the product ((H3 O)+ (MetOH)3 Cl- ) can be only few years. Such reaction can take place in Sagittarius B2, where HCl, H2 O and MetOH have already been identified and it can be an important source for the formation of more complex prebiotic structures.

Photochemistry of Monohydrated Chloromethane: Formation of Free and Hydrated Cl- and CH3+ Ions from a Solvent-Shared Semi-Ion-Pair.

The effect of water molecule on the excited states of CH3Cl(H2O), as compared to those of the isolated chloromethane, has been studied at the multireference configuration interaction with singles and doubles (MR-CISD), including extensivity corrections. Eight new Rydberg states are due to the water molecule but the common states of both systems are not severely altered. Potential energy curves of 23 singlet states along the C-Cl coordinate have also been computed at the MR-CISD level. The dissociation energy of the C-Cl bond decreases from ∼0.4 to 0.5 eV due to the water molecule. As for CH3Cl (de Medeiros, V. C., J. Am. Chem. Soc. 2016, 138, 272-280), a stable ion-pair has also been characterized. However, for CH3Cl(H2O), this ion-pair is better described as a solvent-shared semi-ion-pair, CH3+δ(H2O)Cl-δ. This species is connected with three ionic dissociation channels, with two being due to the water molecule. The presence of these new ionic channels, particularly the lowest energy one, [H3C-O]+ + Cl-, raises a very important question of atmospheric relevance: can the interaction of chloroalkanes with water decrease its deleterious effect on the ozone layer? Several potentially new competing dissociation channels are also studied. The latter results can help to set up the most important states to be included in nonadiabatic dynamic calculations to study how the yields of the ionic channels change due to the water molecule.

Can a gas phase contact ion pair containing a hydrocarbon carbocation be formed in the ground state?

So far, no conclusive evidence of a ground-state contact ion-pair containing a hydrocarbon carbocation has been given in the gas phase. Due to the very high stability of the 1,2:4,5-dibenzotropylium (or dibenzo[a,d]tropylium) carbocation, we suggest (supported by DFT and MP2 calculations) the formation of a contact ion pair between this carbocation and chloride, occurring during the reaction between 1,2:4,5-dibenzotropyl (also named dibenzo[a,d]tropyl or dibenzo[a,d]cycloheptenyl) radical and chlorine atom at very low temperatures, through the harpoon mechanism. This is the first modeling study to find computational evidence for the possibility of a gas-phase contact ion pair (containing a hydrocarbon carbocation) formed in the ground state. Identification of this metastable species can be carried out by trapping it in He nanodroplets, along with infrared laser spectroscopy routinely coupled with this technique.

The generality of the GUGA MRCI approach in COLUMBUS for treating complex quantum chemistry.

The core part of the program system COLUMBUS allows highly efficient calculations using variational multireference (MR) methods in the framework of configuration interaction with single and double excitations (MR-CISD) and averaged quadratic coupled-cluster calculations (MR-AQCC), based on uncontracted sets of configurations and the graphical unitary group approach (GUGA). The availability of analytic MR-CISD and MR-AQCC energy gradients and analytic nonadiabatic couplings for MR-CISD enables exciting applications including, e.g., investigations of π-conjugated biradicaloid compounds, calculations of multitudes of excited states, development of diabatization procedures, and furnishing the electronic structure information for on-the-fly surface nonadiabatic dynamics. With fully variational uncontracted spin-orbit MRCI, COLUMBUS provides a unique possibility of performing high-level calculations on compounds containing heavy atoms up to lanthanides and actinides. Crucial for carrying out all of these calculations effectively is the availability of an efficient parallel code for the CI step. Configuration spaces of several billion in size now can be treated quite routinely on standard parallel computer clusters. Emerging developments in COLUMBUS, including the all configuration mean energy multiconfiguration self-consistent field method and the graphically contracted function method, promise to allow practically unlimited configuration space dimensions. Spin density based on the GUGA approach, analytic spin-orbit energy gradients, possibilities for local electron correlation MR calculations, development of general interfaces for nonadiabatic dynamics, and MRCI linear vibronic coupling models conclude this overview.

Photochemistry of CF3Cl: Quenching of Charged Fragments Is Caused by Nonadiabatic Effects.

For the first time, high-level multireference electronic structure calculations have been performed to study the photochemistry of CF3Cl, allowing a comprehensive interpretation and assignment of experimental data concerning fluorescence, ion-pair formation, and generation of CF3 fragments in several electronic states. All studied dissociation channels correlate either with Cl or Cl- in the ground state. On the other hand, a CF3 fragment can be generated either in the ground or excited state. A rationalization for the nonadiabatic relaxation of CF3Cl, including the formation of an (n4s) stable state and internal conversion at multiple-state intersections, has been provided. Our results explain the anomalous quenching of a charged fragment after low-energy excitation, a fact experimentally observed by separate groups. We show that the CF3+···Cl- ion pair undergoes an internal conversion to the ground state, producing neutral CF3 and Cl fragments. The results also allow understanding as to why CF3Cl is usually a nonemitting species and how UV emission could be induced.

Spin-Forbidden Branching in the Mechanism of the Intrinsic Haber-Weiss Reaction.

The mechanism of the O2⋅- and H2O2 reaction (Haber-Weiss) under solvent-free conditions has been characterized at the DFT and CCSD(T) level of theory to account for the ease of this reaction in the gas phase and the formation of two different set of products (Blanksby et al., Angew. Chem. Int. Ed. 2007, 46, 4948). The reaction is shown to proceed through an electron-transfer process from the superoxide anion to hydrogen peroxide, along two pathways. While the O3⋅- + H2O products are formed from a spin-allowed reaction (on the doublet surface), the preferred products, O⋅-(H2O)+3O2, are formed through a spin-forbidden reaction as a result of a favorable crossing point between the doublet and quartet surface. Plausible reasons for the preference toward the latter set are given in terms of the characteristics of the minimum energy crossing point (MECP) and the stability of an intermediate formed (after the MECP) in the quartet surface. These unique results show that these two pathways are associated with a bifurcation, yielding spin-dependent products.

Challenges encountered during development of Mn porphyrin-based, potent redox-active drug and superoxide dismutase mimic, MnTnBuOE-2-PyP5+, and its alkoxyalkyl analogues.

We disclose here the studies that preceded and guided the preparation of the metal-based, redox-active therapeutic Mn(III) meso-tetrakis(N-n-butoxyethylpyridyl)porphyrin, MnTnBuOE-2-PyP5+ (BMX-001), which is currently in Phase I/II Clinical Trials at Duke University (USA) as a radioprotector of normal tissues in cancer patients. N-substituted pyridylporphyrins are ligands for Mn(III) complexes that are among the most potent superoxide dismutase mimics thus far synthesized. To advance their design, thereby improving their physical and chemical properties and bioavailability/toxicity profiles, we undertook a systematic study on placing oxygen atoms into N-alkylpyridyl chains via alkoxyalkylation reaction. For the first time we show here the unforeseen structural rearrangement that happens during the alkoxyalkylation reaction by the corresponding tosylates. Comprehensive experimental and computational approaches were employed to solve the rearrangement mechanism involved in quaternization of pyridyl nitrogens, which, instead of a single product, led to a variety of mixed N-alkoxyalkylated and N-alkylated pyridylporphyrins. The rearrangement mechanism involves the formation of an intermediate alkyl oxonium cation in a chain-length-dependent manner, which subsequently drives differential kinetics and thermodynamics of competing N-alkoxyalkylation versus in situ N-alkylation. The use of alkoxyalkyl tosylates, of different length of alkyl fragments adjacent to oxygen atom, allowed us to identify the set of alkyl fragments that would result in the synthesis of a single compound of high purity and excellent therapeutic potential.

Revisiting the concept of the (a)synchronicity of Diels-Alder reactions based on the dynamics of quasiclassical trajectories.

A number of model Diels-Alder (D-A) cycloaddition reactions (H2C=CH2 + cyclopentadiene and H2C=CHX + 1,3-butadiene, with X = H, F, CH3, OH, CN, NH2, and NO) were studied by static (transition state - TS and IRC) and dynamics (quasiclassical trajectories) approaches to establish the (a)synchronous character of the concerted mechanism. The use of static criteria, such as the asymmetry of the TS geometry, for classifying and quantifying the (a)synchronicity of the concerted D-A reaction mechanism is shown to be severely limited and to provide contradictory results and conclusions when compared to the dynamics approach. The time elapsed between the events is shown to be a more reliable and unbiased criterion and all the studied D-A reactions, except for the case of H2C=CHNO, are classified as synchronous, despite the gradual and quite distinct degrees of (a)symmetry of the TS structures.

Photochemistry of CH3Cl: Dissociation and CH···Cl Hydrogen Bond Formation.

State-of-the-art electronic structure calculations (MR-CISD) are used to map five different dissociation channels of CH3Cl along the C-Cl coordinate: (i) CH3(X̃(2)A2″) + Cl((2)P), (ii) CH3(3s(2)A1') + Cl((2)P), (iii) CH3(+)((1)A1') + Cl(-)((1)S), (iv) CH3(3p(2)E') + Cl((2)P), and (v) CH3(3p(2)A2″) + Cl((2)P). By the first time these latter four dissociation channels, accessible upon VUV absorption, are described. The corresponding dissociation limits, obtained at the MR-CISD+Q level, are 3.70, 9.50, 10.08, 10.76, and 11.01 eV. The first channel can be accessed through nσ* and n3s states, while the second channel can be accessed through n(e)3s, n(e)3p(σ), and σ3s states. The third channel, corresponding to the CH3(+) + Cl(-) ion-pair, is accessed through n(e)3p(e) states. The fourth is accessed through n(e)3p(e), n(e)3p(σ), and σ3p(σ), while the fifth through σ3p(e) and σ(CH)σ* states. The population of the diverse channels is controlled by two geometrical spots, where intersections between multiple states allow a cascade of nonadiabatic events. The ion-pair dissociation occurs through formation of CH3(+)···Cl(-)and H2CH(+)···Cl(-) intermediate complexes bound by 3.69 and 4.65 eV. The enhanced stability of the H2CH(+)···Cl(-) complex is due to a CH···Cl hydrogen bond. A time-resolved spectroscopic setup is proposed to detect those complexes.

Assessment of density-functionals for describing the X(-) + CH3ONO2 gas-phase reactions with X = F, OH, CH2CN.

The energetics of the ECO2, SN2@C and SN2@N channels of X(-) + CH3ONO2 (X = F, OH, CH2CN) gas-phase reactions were computed using the CCSD(T)/CBS method. This benchmark extends a previous study with X = OH [M. A. F. de Souza et al., J. Am. Chem. Soc., 2012, 134, 19004] and was used to ascertain the accuracy and robustness of nineteen density-functionals for describing these potential energy profiles (PEP) as well as the kinetic product distributions obtained from RRKM calculations. Assessments were based on the mean unsigned error (MUE), the mean signed error (MSE), the #best : #worst (BW) criterion and the statistical confidence interval (CI) for the MSE. In general, double-hybrid (DH) functionals perform better than the range-separated ones, and both are better than the global-hybrid functionals. Based on the MUE and CI criteria the B2GPPLYP, B2PLYP, M08-SO, BMK, ωB97X-D, CAM-B3LYP, M06, M08-HX, ωB97X and B97-K functionals show the best performance in the description of these PEPs. Within this set, the B2GPPLYP functional is the most accurate and robust. The RRKM results indicate that the DHs are the best for describing the selectivities of these reactions. Compared to CCSD(T), the B2PLYP method has a relative error of only ca. 1% for the selectivity and the accuracy to provide the correct conclusion concerning the nonstatistical behavior of these reactions.

Dynamic effects dictate the mechanism and selectivity of dehydration-rearrangement reactions of protonated alcohols [Me2 (R)CCH(OH2 )Me](+) (R=Me, Et, iPr) in the gas phase.

The gas-phase dehydration-rearrangement (DR) reactions of protonated alcohols [Me2 (R)CCH(OH2 )Me](+) [R=Me (ME), Et (ET), and iPr (I-PR)] were studied by using static approaches (intrinsic reaction coordinate (IRC), Rice-Ramsperger-Kassel-Marcus theory) and dynamics (quasiclassical trajectory) simulations at the B3LYP/6-31G(d) level of theory. The concerted mechanism involves simultaneous water dissociation and alkyl migration, whereas in the stepwise reaction pathway the dehydration step leads to a secondary carbocation intermediate followed by alkyl migration. Internal rotation (IR) can change the relative position of the migrating alkyl group and the leaving group (water), so distinct products may be obtained: [Me(R)CCH(Me)Me⋅⋅⋅OH2 ](+) and [Me(Me)CCH(R)Me⋅⋅⋅OH2 ](+) . The static approach predicts that these reactions are concerted, with the selectivity towards these different products determined by the proportion of the conformers of the initial protonated alcohols. These selectivities are explained by the DR processes being much faster than IR. These results are in direct contradiction with the dynamics simulations, which indicate a predominantly stepwise mechanism and selectivities that depend on the alkyl groups and dynamics effects. Indeed, despite the lifetimes of the secondary carbocations being short (<0.5 ps), IR can take place and thus provide a rich selectivity. These different selectivities, particularly for ET and I-PR, are amenable to experimental observation and provide evidence for the minor role played by potential-energy surface and the relevance of the dynamics effects (non-IRC pathways, IR) in determining the reaction mechanisms and product distribution (selectivity).

Ab initio and DFT conformational study on N-nitrosodiethylamine, (C2H5)2N-N=O.

Ab initio (MP2) and DFT (B3LYP) calculations, using the cc-pVTZ and aug-cc-pVTZ basis sets, have been performed to characterize some stationary points on the ground state potential energy surface of the title molecules. Several properties as, for instance, relative energies, the barriers for NO rotation around the NN bond, NBO charges on O and amino N atoms, as well as the dipole moments, have been calculated and analyzed in the light of the structures found. Both computational levels here employed yield three minima, in which the C(2)NNO frame is 'planar' or 'quasi-planar'. Important correlations between NBO charges and geometric parameters, as well as between some structural features and dipole moments are also discussed. A total of 17 structures have been found for the (C(2)H(5))(2)N-N=O molecule. Two ranges of values have been obtained for the dipole moment, with the largest values occurring for the structures in which the nitrogen lone pair is parallel to the NO group π system. For instance, these two ranges are from ~4.1 to 4.5 D, and from ~1.6 to 2.1 D, at the MP2/cc-pVTZ level. These ranges are consistent with a larger and a smaller contribution of a dipolar resonance structure, respectively. As the method or basis set changes the values of the dipole moments change by at most ~0.23 D.

CASSCF and multireference CI with singles and doubles study of low-lying valence and Rydberg states of 2H-tetrazole.

Complete active space self-consistent field (CASSCF) and multireference CI with singles and doubles (MR-CISD) calculations [including extensivity corrections, at MR-CISD+Q and multireference averaged quadratic coupled cluster (MR-AQCC) levels] have been performed to characterize the low-lying valence and the Rydberg states of 2H-tetrazole. The highest level results (MR-AQCC/d'-aug'-cc-pVDZ) indicate the following ordering of the valence singlet excited states: S(1) (n-pi*), 6.06 eV; S(2) (n-pi*), 6.55 eV; S(3) (pi-pi*), 6.55 eV. The MR-CISD+Q/d'-aug'-cc-pVDZ results indicate the same ordering, but at slight higher energies: 6.16, 6.68, and 6.69 eV, respectively. According to our MR-CISD+Q/d'-aug'-cc-pVDZ results, the next two states are Rydberg states, at 7.69 eV (pi-3s) and 7.89 eV (n-3s). The calculated energies of these two states, as well as their proximity, are consistent with the conclusion reached by Palmer and Beveridge (Chem Phys 1987, 111, 249) that the first band of the photoelectron spectrum of 2H-tetrazole is likely to be associated to the first two ionizations processes (of pi and N lone pair electrons), at energies close to 11.3 eV.

Dissociation of ground and nsigma* states of CF3Cl using multireference configuration interaction with singles and doubles and with multireference average quadratic coupled cluster extensivity corrections.

Extended complete active space self-consistent field (CASSCF), multireference configuration interaction with singles and doubles (MR-CISD), and multireference average quadratic coupled cluster (MR-AQCC) calculations have been performed on the ground (S(0)) and first excited (nsigma(*),S(1)) states of the CF(3)Cl molecule. Full geometry optimizations have been carried out for S(0) as well as "relaxed" potential energy calculations for both states, along the C-Cl bond distance. Vertical excitation energies (DeltaE(vertical)), dissociation energies (DeltaE(diss)), dissociation enthalpies (DeltaH(diss)), and the oscillator strength (f) have also been computed. Basis set effects, basis set superposition error (BSSE), and spin-orbit and size-extensivity corrections have also been considered. The general agreement between theoretical and available experimental results is very good. The best results for the equilibrium geometrical parameters of S(0) (at MR-AQCCaug-cc-pVTZ+d level) are 1.762 and 1.323 A, for the C-Cl and C-F bond distances, respectively, while the corresponding experimental values are 1.751 and 1.328 A. The [angle](ClCF) and [angle](FCF) bond angles are in excellent agreement with the corresponding experimental values (110.3 degrees and 108.6 degrees ). The best calculated values for DeltaE(vertical), DeltaH(diss), and f are 7.63 eV [at the MR-AQCCaug-cc-pV(T+d)Z level], 3.59 eV[MR-AQCCaug-cc-pV(T+d)Z level+spin-orbit and BSSE corrections], and 2.74x10(-3) (MR-CISD/cc-pVTZ), in comparison with the corresponding experimental values of 7.7+/-0.1 eV, 3.68 eV, and 3.12 x 10(-3)+/-2.50 x 10(-4). The results concerning the potential energy curves for S(0) and S(1) show a tendency toward the nonoccurrence of crossing between these two states (in the intermediate region along the C-Cl coordinate), as the basis set size increases. Such tendency is accompanied by a decreasing well depth for the S(1) state. Dynamic electronic correlation (especially at the MR-AQCC level) is also an important factor toward an absence of crossing along the C-Cl coordinate. Further investigations of a possible crossing using gradient driven techniques (at CASSCF and MR-CISD levels) seem to confirm its absence.

An ab initio study of the C2H2-HF, C2H(CH3)-HF and C2(CH3)2-HF hydrogen-bonded complexes.

MP2/6-31++G** and B3LYP/6-31++G** ab initio molecular orbital calculations have been performed in order to obtain molecular geometries, binding energies and vibrational properties of the C2H2-HF, C2H(CH3)-HF and C2(CH3)2-HF H-bonded complexes. As expected, the more pronounced effects on the structural properties of the isolated molecules due to complexation was verified for the C[triple bond]C and H-F bond lengths, which are directly involved in the H-bond formation. These bond distances increased after complexation. BSSE uncorrected B3LYP binding energies are always lower than the corresponding MP2 values. However, the opposite trend has been verified after BSSE correction by the counterpoise method since it is much lower at B3LYP than at MP2 level. The binding energies for these complexes as well as for the HF acid submolecule modes (the HF stretching and vibrational frequency modes) showed an increasing hydrogen-bonding strength with increasing methyl substitution. The splitting in the HF in-plane and out-of-plane bending modes reflects the anisotropy in the hydrogen-bonding interaction with the pi system of the C[triple bond]C bond. The H-F stretching frequency is shifted downward after complexation and it increases with the methyl substitution. The IR intensities of the HF acid submolecule fundamentals are adequately interpreted through the atomic polar tensor of the hydrogen atom using the charge-charge flux-overlap model. The skeletal stretching modes of the Alkyne submolecule are decreased in the complex. The new vibrational modes arising from complexation show several interesting features.

A theoretical study of hydrogen complexes of the XH-pi type between propyne and HF, HCL or HCN.

The present manuscript reports a systematic investigation of the basis set dependence of some properties of hydrogen-bonded (pi type) complexes formed by propyne and a HX molecule, where X=F, Cl and CN. The calculations have been performed at Hartree-Fock, MP2 and B3LYP levels. Geometries, H-bond energies and vibrational have been considered. The more pronounced effects on the structural parameters of the isolated molecules, as a result of complexation, are verified on RCtriple bondC and HX bond lengths. As compared to double-zeta (6-31G**), triple-zeta (6-311G**) basis set leads to an increase of RCtriple bondC bond distance, at all three computational levels. In the case where diffuse functions are added to both hydrogen and 'heavy' atoms, the effect is more pronounced. The propyne-HX structural parameters are quite similar to the corresponding parameters of acetylene-HX complexes, at all levels. The largest difference is obtained for hydrogen bond distance, RH, with a smaller value for propyne-HX complex, indicating a stronger bond. Concerning the electronic properties, the results yield the following ordering for H-bond energies, DeltaE: propynecdots, three dots, centeredHF>propynecdots, three dots, centeredHCl>propynecdots, three dots, centeredHCN. It is also important to point out that the inclusion of BSSE and zero-point energies (ZPE) corrections cause significant changes on DeltaE. The smaller effect of ZPE is obtained for propynecdots, three dots, centeredHCN at HF/6-311++G** level, while the greatest difference is obtained at MP2/6-31G** level for propynecdots, three dots, centeredHF system. Concerning the IR vibrational it was obtained that larger shift can be associated with stronger hydrogen bonds. The more pronounced effect on the normal modes of the isolated molecule after the complexation is obtained for HX stretching frequency, which is shifted downward.